Hidden Orbital Order in URU2S12 
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When matter is cooled from high temperatures, collective instabilities de- 
velop amongst its constituent particles that lead to new kinds of order.S An 
anomaly in the specific heat is a classic signature of this phenomenon. Usu- 
ally the associated order is easily identified, but sometimes its nature remains 
elusive. The heavy fermion metal URU2S12 is one such example, where the or- 
der responsible for the sharp specific heat anomaly at T = 17K has remained 
unidentified despite more than seventeen years of effort.i In URU2S12, the coex- 
istence of large electron-electron repulsion and antiferromagnetic fluctuations in 
URU2S12 leads to an almost incompressible heavy electron fluid, where anisotrop- 
ically paired quasiparticle states are energetically favored.! In this paper we use 
these insights to develop a detailed proposal for the hidden order in URU2S12. 
We show that incommensurate orbital antiferromagnetism, associated with cir- 
culating currents between the uranium ions, can account for the local fields 
and entropy loss observed at the 17K transition; furthermore we make detailed 
predictions for neutron scattering measurements. 

The intermetallic compound URU2S12 contains a dense lattice of local moments where 
quantum fluctuations barely prevent spin ordering; the residual antiferromagnetic couplings 
between the strongly repulsive electrons are large and can drive new types of collective 
instabilities.! At To = 17.5K, URU2S12 undergoes a second-order phase transition char- 



acterized by sharp discontinuities in bulk properties including specific heat,! lineaiH and 
nonlinear!!! susceptibilities, thermal expansiorJll and resistivity.! The accompanying gap in 
the magnetic excitation spectrumi suggests the formation of an itinerant spin density wave 
at this temperature; however the size of the observed staggered moment^ (mo = 0.03/i£) 
cannot account for the entropy loss at this transition. The distinction between the primary 
hidden order parameter and the secondary magnetic order parameter is clarified by high- 
field measurement s0'0 which indicate that the bulk anomalies survive up to an applied 
field strength of 40T, whereas the magnetically ordered moment is destroyed by fields less 
than half this magnitude (15T).i! Furthermore the size of the small ordered moment grows 
linearly with pressure^ while the gap associated with the hidden order is relatively robust 
over this pressure range.B 

There have been many theoretical proposals for the primary order parameter in 
U 7?-u 2 S22-0 Two recent NMR studies have provided new insights on this outstanding prob- 
lem with crucial consequences for the nature of the hidden order. It has been widely assumed 
that the spin antiferromagnetism and the hidden order are coupled and homogeneous. How- 
ever a recent NMR study of URu 2 Si 2 under pressure^ indicates that for T < T there 
exist distinct antiferromagnetic and paramagnetic regions, implying that the magnetic and 
the hidden orders are phase-separated. This conclusion, supported by /jSR data,@ implies 
that the hidden order phase contains no spin order. The observed growth of the staggered 
spin moment with applied pressure^ is then simply a volume fraction effect which develops 
separately from the hidden order via a first order transition. At ambient pressure roughly 
a tenth or less of the system0'@ is magnetic with m sp i n = 0.3/ie. The key implication of 
these measurements for theory is that the magnetic and the hidden order parameters are 
independent .0 

In a parallel study, NMR measurements at ambient pressure^ on URu2Si2 indicate 
that at T < To the central (non-split) silicon NMR line-width develops a field-independent, 
isotropic component whose temperature-dependent magnitude is proportional to that of the 
hidden order parameter. These results imply an isotropic field distribution at the silicon 



sites whose root-mean square value is proportional to the hidden order (ip) 



(1) 



and is ~ 10 Gauss at T = 0. This field magnitude is too small to be explained by the 



bond length (a = 4 x 10 -8 cm). Furthermore this moment is aligned along the c-axis, and 
thus cannot account for the isotropic nature of the local field distribution detected by NMR. 
These measurements indicate that as the hidden order develops, an isotropically distributed 
static magnetic field develops at each silicon site. This is unambiguous evidence that the 
hidden order parameter breaks time-reversal invariance. 

Guided by these recent experiments, we now discuss our proposal for the nature of the 
hidden order parameter. The magnetic fields at the silicon nuclei have two possible origins^ 
the conduction electron-spin interaction and the orbital shift due to current densities. In 
URu 2 Si2, the electron fluid exhibits a strong Ising anisotropy along the c-axis, as measured 
by the Knight shift ;0 hence electron-spin coupling cannot be responsible for the observed 
isotropic fields. It would thus appear that these local fields are induced by currents that 
develop inside the crystal as the hidden order develops, and accordingly we attribute the 
observed isotropic linewidth to the orbital shift. 

In this paper we propose that URU2S12 becomes an incommensurate orbital antifer- 
romagnet (OAF) at T = T with charge currents circulating between the uranium ions. 
Similar states have been studied extensively in the context of the two-dimensional Hubbard 
model,iHHl particularly in connection with staggered flux phases.EHH More recently com- 
mensurate current density wave order has been proposed to explain the spin-gap phase in 
the underdoped cuprate superconductors.^^ The planar tetragonal structure of URu 2 Si2 
lends itself naturally to an anisotropic charge instability of this type. Here we show that 
incommensurate orbital antiferromagnetic order in \JRu2Si2 can quantitatively account for 
the existing specific heat and NMR data for T < To. In order to test this proposal, we 
make specific predictions for neutron scattering. Calculation of the structure factor re- 
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quires knowledge of the fields throughout the full volume in real-space. Since the NMR 
only yields this information at discrete points we need some additional input to proceed. 
In the orbital antiferromagnet, the spatial dependence of the fields throughout the sample 
is determined by the requirement that the field distribution at the silicon sites is isotropic. 
Using this approach, we are able to link quantitatively the fields observed by NMR to the 
large specific heat anomaly that develops at To. We also use the spatial field distribution as- 
sociated with the incommensurate orbital antiferromagnet to predict the position, intensity 
and form-factor associated with neutron scattering peaks on a surface of constant anisotropy 
in momentum space centered on the origin. 

We begin by estimating the local fields at the silicon sites due to orbital currents circu- 
lating around the square uranium plaquettes in the a — b plane. On dimensional grounds, 
the current along the U — U bond is given by / ~ ^ where A is the gap associated with 
the formation of hidden order. Then the field induced at a height a above a plaquette 
is B ~ f^M {jit) = 11 Gauss, in good agreement with the local field strength detected 
in NMR measurements; here&i~0 we have used A = 11QK. The resulting orbital mo- 
ment, m or b ~ 0.02/is, is comparable to the effective spin moment at ambient pressure 
(m = m spin [P amb ] = [.10}m spin = 0.03/xb). We emphasize that an orbital moment produces 
a local field an order of magnitude less than that associated with a spin moment of the same 
value; the low field strengths observed at the silicon sites are quantitatively consistent with 
our proposal that they originate from charge currents. 

An orbital moment, m or b = 0.02fi B , can also account for the sizeable entropy loss at the 
transition. In a metal the change in entropy is given by AS = A7„T , where A7„ is the 
change in the linear specific heat coefficient resulting from the gapping of the Fermi surface. 
A~f n is inversely proportional to the Fermi energy ep of the gapped Fermi surface. Thus 
in general, the change in entropy per unit cell is given by AS ~ k B ^^. Exploiting the 
mean-field nature of this transition (2 A oc To), we find that 

ASoafm ~ k B = 0.02k B (2) 
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where fi B = (j^J a 2 ep, is the saturated orbital moment reached when the hidden order 
gap A is approximately ep. The ratio \ibI\i*q = (^f) 2 (ffr)> where ao is the Bohr radius 
and €h = e 2 /(2ao) is the energy of a hydrogen atom. For URu2Si2, the very large size 
of \t^\ ~ 10 3 produced by the large mass renormalization of the heavy electrons actually 
offsets the small ratio ^ ~ 10 _1 , so that ^b/ ~ 10, and ASoafm ~ 0.2ks = 0.3ks In 2 is a 
number in good agreement with experiment .@ The critical field for suppressing the associated 

jjorb 

thermodynamic anomalies is distinct from its spin counterpart; the ratio ? pin ~ ~ 10 
is qualitatively consistent with the observed critical field of approximately 40T associated 
with the destruction of the hidden orderSHiil. From this simple discussion, we see that 
the sizable entropy loss associated with the development of the orbital antiferromagnetic 
state is a direct consequence of the strong renormalization of the electron mass in URU2S12 
^ ef^' ^ e absence of such a large effective mass (^- ~ 3) in the cuprates could explain 
why analogous thermodynamic anomalies are difficult to observe there. 

Next we test whether the proposed orbital antiferromagnetism will yield the observed 
isotropic local fields. In order to do so, we allow the circulating current around a plaquette 
(cf. Fig. 1) centered at site X to develop a modulated magnetization M(X) = ipe lCi x . 
The current along a bond is then the difference of the circulating currents along its adjacent 
plaquettes. The field at a silicon site can be computed using Ampere's law, where the 
relevant vector potential is 



1 _ f*f , , Kx 



A(x) = :E ^ dx' j- T 3 \^r (3) 



cyj^p |x-( Xj + x'). 

where xi- ' are the endpoints of the bond at site Xj. 

The silicon atoms in URU2S12 are located at low-symmetry sites, so that the fields do not 
cancel; they reside above and below the centers of the uranium plaquettes. Microscopically, 
wave-vector selection is most likely due to details of the Fermi surface, however we can 
obtain a good idea of the likely modulation Q vector from the isotropic nature of the field 
distribution at the silicon sites. For example, in the commensurate case Q = (1/2, 1/2, 1), 
the fields on the silicon sites are only along the c-axis and thus the field distribution would 



be highly anisotropic. Consider a wave-vector (q, q, 1) (c.f. Fig. 1). In this case, the currents 
are modulated within a plane with a wavelength 2n/q, but staggered between planes. This 
then produces a circulating field distribution where the component of the field parallel to 
the (0,0, 1) planes becomes larger and larger as q is reduced. To obtain an isotropic field 
distribution, q needs to be reduced to a point where the horizontal and vertical components 
of the field are comparable. A detailed calculation based on the above model, shows that 
the condition of perfect isotropy yields a circle of wave- vectors (Fig. 2.) centered around 
Q = (0, 0, 1) with a radius q ss 0.22. Relaxation of this constraint results in an annulus of 
possible Q vectors as shown in Fig. 2. 

Our proposal of incommensurate current ordering in URu 2 Si2 can be tested by experi- 
ment. In particular we can Fourier transform the real-space magnetic fields to calculate the 
neutron scattering cross-section 

£-(£)W-*tt>- w 

Here \B(q) | 2 is the structure factor of the magnetic fields produced by the orbital currents and 
S(q) = \B(q)\ 2 /(4:7iii B ) 2 is the structure factor of the orbital magnetic moments, measured 
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in units of electron Bohr magnetons (/ie). The factor r = ^ e c2 , m e is the electron mass 
and (7jv is the neutron gyromagnetic ratio. Using the vector potential from eq. (3), we have 
calculated the magnetic field distribution for the incommensurate orbital antiferromagnet 
described above, and find that its Fourier transform is given by 

B(q) = ^Nla* £ 5 q , Q+G Jof^U^l (1 + (-!)(-+-+-)) • f^M] x q, (5) 
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where jo(x) = From this expression, we can determine the intensities and the form 

factors associated with the diffraction. We find that there exists a set of dominant peaks 
associated with a constant anisotropy locus of wavevectors (Fig. 2.) in the first Brillouin 
zone (BZ) where S(q) = 0.18 (^) , indicating that roughly a fifth of the total integrated 
weight of S(q) (TIW) lies here. Using the sum rule that relates the TIW to the square 
of the moment, and the fact that at ambient pressure only a tenth of the sample is (spin) 
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magnetic, we find that the scattering peaks due to orbital ordering in the first BZ should 
have 1/50 the intensity of the analogous spin magnetic peaks at ambient pressure. We have 
also calculated that these peaks will have a form-factor (cf. inset Fig. 2) that scales with 
wavevector as g~ 4 , where this power-law decay is signatory of an extended scattering source. 

We end with remarks about the microscopic formation of these charge currents. At low 
temperatures heavy electron materials form almost incompressible Fermi liquids. The large 
Coulomb repulsion between the Landau quasiparticles strongly suppresses on-site charge 
fluctuations, demanding nodes in the particle-hole wavefunction. The resulting anisotropi- 
cally paired states are also favored by the residual antiferromagnetic interactions that persist 
in the heavy electron fluid. Indeed we believe that the same d-wave pairing that drives the 
superconducting transition in URu 2 Si2 at T = 1.5K also plays an important role in the 
formation of the orbital antiferromagnetism. We have found that the development of such 
anisotropic charge-density wave pairing occurs naturally in a simple model of URU2S12 with 
nearest-neighbor antiferromagnetic interactions, 

H I = ^Aq)S(q)-S(-q). (6) 

q 

Here S(q) is the Fourier transform of the magnetization at wavevector q and J(q) = 
2J(cosq x + cosq y ) for a nearest- neigb our interaction between adjacent U atoms in the basal 
plane. Expanding this interaction in terms of quasiparticle operators, we find that it can be 
re-written as a sum of attractive interactions in four independent anisotropic charge density 
wave channels: 

Hi = -J E (t£)*t£, p kl (Q)PkA-Q)- (?) 

r=i,4,fci,fe2,Q 

Here pt{Q) = J2a=±i/2 c \-Q/2a c k+Q/2a is the charge density operator at wavevector Q, ex- 
pressed in terms of creation and annihilation operators of the heavy quasiparticles. The 
form factors are 'y 1,2 (k) = cos(k x ) ± cos(k y ),^ 3 ' 4 (k) = i{sm(k x ) ±sm(k y )). Of these four 
possibilities, channels 1 and 3 do not have nodes, and are therefore suppressed by local 
Coulomb interactions. In the remaining two channels, only T = 4 breaks time reversal 
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symmetry, giving rise to incommensurate orbital currents. One of the questions for further 
study concerns the excitations of this mean-field state. In particular, an incommensurate 
OAF is expected to exhibit a gapless phason mode associated with slow fluctuations of its 
wavevector. URu2Si2 is known to develop a dispersing singlet excitation^ at T = To. This 
mode was previously attributed to spin antiferromagnetism, now known to be absent from 
the hidden order phase.0 We plan to study the possible identification of this propagating 
mode with the phason of an incommensurate orbital antiferromagnet. 

In summary, we have discussed the theoretical implications of two recent NMR exper- 
iments on the hidden order in URu2Si2- Pressure-dependent measurements indicate that 
it is completely independent of the spin magnetism in this material. We argue that the 
development of isotropically distributed magnetic fields at the silicon sites at T = T implies 
that the hidden order parameter breaks time-reversal symmetry. The size and the isotropy 
of these observed local fields lead us to propose that URU2S12 becomes an incommensurate 
orbital antiferromagnet at T < T . The heavy electron mass reduces the saturation value 
of the orbital moment, accounting for the sizable entropy loss at the transition and the 
scale of the associated critical field. We calculate the positions, intensities and form-factor 
associated with the resulting neutron scattering peaks so that this proposal can be tested 
by experiment. 
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Figure Captions 

Fig. 1. Magnetic field distribution associated with incommensurate orbital currents in the 
(0, 0, 1) plane, (a) Schematic illustration of incommensurate orbital currents, showing 
resulting magnetic fields above and below the (0, 0, 1) plane, (b) Side view showing 
how proposed field distribution is isotropic at silicon sites with a staggered current 
distribution between layers corresponding to Q — (0.16,0.16, 1). 

Fig. 2. An isotropic magnetic field distribution at the silicon sites can be produced by many 
different wave vectors Q for orbital order, all of which give qualitatively similar neutron 
scattering patterns, (a) Contour plot showing locus of constant anisotropy around 
a ring of radius Q± = 0.22 in the vicinity of Q = (0,0,1). (b) Predicted elastic 
neutron scattering intensity, where S(q)% = S(q)%/ (^~) ~ \]\B(q)\ 2 gives a measure 
of the Fourier spectrum of magnetic fields B(q), plotted for the representative case 
Q = (0.16,0.16,1). 
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Fig. 2 
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